Code Interpreter for Bioinformatics: Are We There Yet?

The Code Interpreter feature in ChatGPT has the potential to democratize data analysis for non-specialists. As bioinformaticians, we are impressed by its performance in data manipulation and visualization. However, bioinformatics tasks often require execution of third-party packages, access to annotation knowledgebase, and handling large datasets. Code Interpreter's exclusive support for Python, no installation option for additional packages, inability to utilize external resources, and limited storage capacity could pose obstacles to its wide adoption in bioinformatics applications. To address these limitations, we advocated for the necessity of locally deployable, API-based systems for chatbot-aided bioinformatics applications.

Subfamily-specific differential contribution of individual monomers and the tether sequence to mouse L1 promoter activity.

BACKGROUND: The internal promoter in L1 5'UTR is critical for autonomous L1 transcription and initiating retrotransposition. Unlike the human genome, which features one contemporarily active subfamily, four subfamilies (A_I, Gf_I and Tf_I/II) have been amplifying in the mouse genome in the last one million years. Moreover, mouse L1 5'UTRs are organized into tandem repeats called monomers, which are separated from ORF1 by a tether domain. In this study, we aim to compare promoter activities across young mouse L1 subfamilies and investigate the contribution of individual monomers and the tether sequence. RESULTS: We observed an inverse relationship between subfamily age and the average number of monomers among evolutionarily young mouse L1 subfamilies. The youngest subgroup (A_I and Tf_I/II) on average carry 3-4 monomers in the 5'UTR. Using a single-vector dual-luciferase reporter assay, we compared promoter activities across six L1 subfamilies (A_I/II, Gf_I and Tf_I/II/III) and established their antisense promoter activities in a mouse embryonic fibroblast cell line and a mouse embryonal carcinoma cell line. Using consensus promoter sequences for three subfamilies (A_I, Gf_I and Tf_I), we dissected the differential roles of individual monomers and the tether domain in L1 promoter activity. We validated that, across multiple subfamilies, the second monomer consistently enhances the overall promoter activity. For individual promoter components, monomer 2 is consistently more active than the corresponding monomer 1 and/or the tether for each subfamily. Importantly, we revealed intricate interactions between monomer 2, monomer 1 and tether domains in a subfamily-specific manner. Furthermore, using three-monomer 5'UTRs, we established a complex nonlinear relationship between the length of the outmost monomer and the overall promoter activity. CONCLUSIONS: The laboratory mouse is an important mammalian model system for human diseases as well as L1 biology. Our study extends previous findings and represents an important step toward a better understanding of the molecular mechanism controlling mouse L1 transcription as well as L1's impact on development and disease.

Influence of conceptus presence and preovulatory estradiol exposure on uterine gene transcripts and proteins around maternal recognition of pregnancy in beef cattle.

The uterine environment must provide sufficient endocrine conditions and nutrients for pregnancy maintenance and conceptus survival. The objective of this study was to determine the effects of preovulatory estradiol and conceptus presence on uterine transcripts and uterine luminal fluid (ULF) proteins. Beef cows/heifers were synchronized and artificially inseminated (d 0). Uteri were flushed (d 16); conceptuses and endometrial biopsies were collected. Total cellular RNA was extracted from endometrium for RNA sequencing and RT-PCR validation. There were two independent ULF pools made for each of the following groups: highE2/conceptus, highE2/noconceptus, lowE2/conceptus, and lowE2/noconceptus that were analyzed using the 2D LC-MS/MS based iTRAQ method. There were 64 differentially expressed genes (DEGs) and 77 differentially expressed proteins (DEPs) in common among the highE2/conceptus vs highE2/noconceptus and lowE2/conceptus vs lowE2/noconceptus groups. In summary, the interaction between preovulatory estradiol and the conceptus induces the expression of genes, proteins, and pathways necessary for pregnancy.

nCov2019: an R package for studying the COVID-19 coronavirus pandemic.

BACKGROUND: The global spreading of the COVID-19 coronavirus is still a serious public health challenge. Although there are a large number of public resources that provide statistics data, tools for retrospective historical data and convenient visualization are still valuable. To provide convenient access to data and visualization on the pandemic we developed an R package, nCov2019 (https://github.com/YuLab-SMU/nCov2019). METHODS: We collect stable and reliable data of COVID-19 cases from multiple authoritative and up-to-date sources, and aggregate the most recent and historical data for each country or even province. Medical progress information, including global vaccine development and therapeutics candidates, were also collected and can be directly accessed in our package. The nCov2019 package provides an R language interfaces and designed functions for data operation and presentation, a set of interfaces to fetch data subset intuitively, visualization methods, and a dashboard with no extra coding requirement for data exploration and interactive analysis. RESULTS: As of January 14, 2021, the global health crisis is still serious. The number of confirmed cases worldwide has reached 91,268,983. Following the USA, India has reached 10 million confirmed cases. Multiple peaks are observed in many countries. Under the efforts of researchers, 51 vaccines and 54 drugs are under development and 14 of these vaccines are already in the pre-clinical phase. DISCUSSION: The nCov2019 package provides detailed statistics data, visualization functions and the Shiny web application, which allows researchers to keep abreast of the latest epidemic spread overview.

Influence of estradiol on bovine trophectoderm and uterine gene transcripts around maternal recognition of pregnancy†.

Embryo survival and pregnancy success is increased among animals that exhibit estrus prior to fixed time-artificial insemination, but there are no differences in conceptus survival to d16. The objective of this study was to determine effects of preovulatory estradiol on uterine transcriptomes, select trophectoderm (TE) transcripts, and uterine luminal fluid proteins. Beef cows/heifers were synchronized, artificially inseminated (d0), and grouped into either high (highE2) or low (lowE2) preovulatory estradiol. Uteri were flushed (d16); conceptuses and endometrial biopsies (n = 29) were collected. RNA sequencing was performed on endometrium. Real-time polymerase chain reaction (RT-PCR) was performed on TE (n = 21) RNA to measure relative abundance of IFNT, PTGS2, TM4SF1, C3, FGFR2, and GAPDH. Uterine fluid was analyzed using 2D Liquid Chromatography with tandem mass spectrometry-based Isobaric tags for relative and absolute quantitation (iTRAQ) method. RT-PCR data were analyzed using the MIXED procedure in SAS. There were no differences in messenger RNA (mRNA) abundances in TE, but there were 432 differentially expressed genes (253 downregulated, 179 upregulated) in highE2/conceptus versus lowE2/conceptus groups. There were also 48 differentially expressed proteins (19 upregulated, 29 downregulated); 6 of these were differentially expressed (FDR < 0.10) at the mRNA level. Similar pathways for mRNA and proteins included: calcium signaling, protein kinase A signaling, and corticotropin-releasing hormone signaling. These differences in uterine function may be preparing the conceptus for improved likelihood of survival after d16 among highE2 animals.

iDEP Web Application for RNA-Seq Data Analysis.

RNA sequencing (RNA-seq) has become a routine method for transcriptomic profiling. We developed a user-friendly web app called iDEP (integrated differential expression and pathway analysis) to help biologists interpret read counts or other types of expression matrices derived from read mapping. With iDEP, users can easily conduct exploratory data analysis, identify differentially expressed genes, and perform pathway analysis. Due to its intuitive user interface and massive annotation database, iDEP is being widely adopted for interactive analysis of RNA-seq data. Using a public dataset on the effect of heat shock on mouse with and without functional Hsf1, we demonstrate how users can prepare data files and conduct in-depth analysis. We also discuss the importance of critical interpretion of results (avoid p-hacking and rationalizing) and validation of significant pathways by using different methods and independent annotation databases.

Hydrotropism in the primary roots of maize.

Recent studies mainly in Arabidopsis have renewed interest and discussion in some of the key issues in hydrotropism of roots, such as the site of water sensing and the involvement of auxin. We examined hydrotropism in maize (Zea mays) primary roots. We determined the site of water sensing along the root using a nonintrusive method. Kinematic analysis was conducted to investigate spatial root elongation during hydrotropic response. Indole-3-acetic acid (IAA) and other hormones were quantified using LC-MS/MS. The transcriptome was analyzed using RNA sequencing. Main results: The very tip of the root is the most sensitive to the hydrostimulant. Hydrotropic bending involves coordinated adjustment of spatial cell elongation and cell flux. IAA redistribution occurred in maize roots, preceding hydrotropic bending. The redistribution is caused by a reduction of IAA content on the side facing a hydrostimulant, resulting in a higher IAA content on the dry side. Transcriptomic analysis of the elongation zone prior to bending identified IAA response and lignin synthesis/wall cross-linking as some of the key processes occurring during the early stages of hydrotropic response. We conclude that maize roots differ from Arabidopsis in the location of hydrostimulant sensing and the involvement of IAA redistribution.

Evaluating the impact of health awareness events on Google search frequency.

Over two hundred health awareness events take place in the United States in order to educate the public about various diseases. It would be informative and instructive for the organizations to know the impact of these events, although such information could be difficult to measure. We investigated whether 46 selected events attract the public attention by increasing the search frequencies of certain keywords. Internet search data from 2004 to 2017 were downloaded from Google Trend (GT). Three statistical methods including Transfer Function Noise modeling, Wilcoxon Rank Sum test, and Binomial inference were conducted. Our study showed that 10 health awareness events resulted in increased search frequencies in the event months, and 28 events did not, with the rest being classified as unclear.

Comparative analysis of methicillin-sensitive and resistant Staphylococcus aureus exposed to emodin based on proteomic profiling.

Emodin has a strong antibacterial activity, including methicillin-resistant Staphylococcus aureus (MRSA). However, the mechanism by which emodin induces growth inhibition against MRSA remains unclear. In this study, the isobaric tags for relative and absolute quantitation (iTRAQ) proteomics approach was used to investigate the modes of action of emodin on a MRSA isolate and methicillin-sensitive S. aureus ATCC29213(MSSA). Proteomic analysis showed that expression levels of 145 and 122 proteins were changed significantly in MRSA and MSSA, respectively, after emodin treatment. Comparative analysis of the functions of differentially expressed proteins between the two strains was performed via bioinformatics tools blast2go and STRING database. Proteins related to pyruvate pathway imbalance induction, protein synthesis inhibition, and DNA synthesis suppression were found in both methicillin-sensitive and resistant strains. Moreover, Interference proteins related to membrane damage mechanism were also observed in MRSA. Our findings indicate that emodin is a potential antibacterial agent targeting MRSA via multiple mechanisms.

ATOH1 Promotes Leptomeningeal Dissemination and Metastasis of Sonic Hedgehog Subgroup Medulloblastomas.

Medulloblastoma arising from the cerebellum is the most common pediatric brain malignancy, with leptomeningeal metastases often present at diagnosis and recurrence associated with poor clinical outcome. In this study, we used mouse medulloblastoma models to explore the relationship of tumor pathophysiology and dysregulated expression of the NOTCH pathway transcription factor ATOH1, which is present in aggressive medulloblastoma subtypes driven by aberrant Sonic Hedgehog/Patched (SHH/PTCH) signaling. In experiments with conditional ATOH1 mouse mutants crossed to Ptch1+/- mice, which develop SHH-driven medulloblastoma, animals with Atoh1 transgene expression developed highly penetrant medulloblastoma at a young age with extensive leptomeningeal disease and metastasis to the spinal cord and brain, resembling xenografts of human SHH medulloblastoma. Metastatic tumors retained abnormal SHH signaling like tumor xenografts. Conversely, ATOH1 expression was detected consistently in recurrent and metastatic SHH medulloblastoma. Chromatin immunoprecipitation sequencing and gene expression profiling identified candidate ATOH1 targets in tumor cells involved in development and tumorigenesis. Among these targets specific to metastatic tumors, there was an enrichment in those implicated in extracellular matrix remodeling activity, cytoskeletal network and interaction with microenvironment, indicating a shift in transcriptomic and epigenomic landscapes during metastasis. Treatment with bone morphogenetic protein or SHH pathway inhibitors decreased tumor cell proliferation and suppressed metastatic tumor growth, respectively. Our work reveals a dynamic ATOH1-driven molecular cascade underlying medulloblastoma metastasis that offers possible therapeutic opportunities. Cancer Res; 77(14); 3766-77. ©2017 AACR.

Root isoflavonoids and hairy root transformation influence key bacterial taxa in the soybean rhizosphere.

Rhizodeposits play a key role in shaping rhizosphere microbial communities. In soybean, isoflavonoids are a key rhizodeposit component that aid in plant defense and enable symbiotic associations with rhizobia. However, it is uncertain if and how they influence rhizosphere microbial communities. Isoflavonoid biosynthesis was silenced via RNA interference of isoflavone synthase in soybean hairy root composite plants. Rhizosphere soil fractions tightly associated with roots were isolated, and PCR amplicons from 16S rRNA gene variable regions V1-V3 and V3-V5 from these fractions were sequenced using 454. The resulting data was resolved using MOTHUR and vegan to identify bacterial taxa and evaluate changes in rhizosphere bacterial communities. The soybean rhizosphere was enriched in Proteobacteria and Bacteroidetes, and had relatively lower levels of Actinobacteria and Acidobacteria compared with bulk soil. Isoflavonoids had a small effect on bacterial community structure, and in particular on the abundance of Xanthomonads and Comamonads. The effect of hairy root transformation on rhizosphere bacterial communities was largely similar to untransformed plant roots with approximately 74% of the bacterial families displaying similar colonization underscoring the suitability of this technique to evaluate the influence of plant roots on rhizosphere bacterial communities. However, hairy root transformation had notable influence on Sphingomonads and Acidobacteria.

Genome-wide survey of interindividual differences of RNA stability in human lymphoblastoid cell lines.

The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.

Regulation of gene expression with thyroid hormone in rats with myocardial infarction.

INTRODUCTION: The expression of hundreds of genes is altered in response to left ventricular (LV) remodeling following large transmural myocardial infarction (MI). Thyroid hormone (TH) improves LV remodeling and cardiac performance after MI. However, the molecular basis is unknown. METHODS: MI was produced by ligation of the left anterior descending coronary artery in female SD rats. Rats were divided into the following groups: (1) Sham MI, (2) MI, and (3) MI+T4 treatment (T4 pellet 3.3 mg, 60 days release, implanted subcutaneously immediately following MI). Four weeks after surgery, total RNA was isolated from LV non-infarcted areas for microarray analysis using the Illumina RatRef-12 Expression BeadChip Platform. RESULTS: Signals were detected in 13,188 genes (out of 22,523), of which the expression of 154 genes were decreased and the expression of 200 genes were increased in MI rats compared with Sham MI rats (false discovery rate (FDR) <0.05). Compared to MI rats, T4 treatment decreased expression of 27 genes and increased expression of 28 genes. In particular, 6 genes down-regulated by MI and 12 genes up-regulated by MI were reversed by T4. Most of the 55 genes altered by T4 treatment are in the category of molecular function under binding (24) and biological processes which includes immune system process (9), multi-organism process (5) and biological regulation (19) nonexclusively. CONCLUSIONS: These results suggest that altered expression of genes for molecular function and biological process may be involved in the beneficial effects of thyroid hormone treatment following MI in rats.

Early gene response of human brain microvascular endothelial cells to Listeria monocytogenes infection.

The gene expression of human brain microvascular endothelial cells (HBMEC) in response to 4 h of infection by Listeria monocytogenes was analyzed. Four hours after infection, the expression of 456 genes of HBMEC had changed (p < 0.05). We noted that many active genes were involved in the formyl-methionyl-leucyl-phenylalanine pathway in infected HBMEC. In the upregulated genes, mRNA levels of interleukin-8 and interleukin-15 in infected cells increased according to microarray and real-time reverse transcription - PCR analyses. Since both cytokines are regarded as potent chemotactic factors, the results suggest that HBMEC are capable of recruiting cells of innate and adaptive immune responses during early L. monocytogenes infection.

Gene expression profiles of estrogen receptor-positive and estrogen receptor-negative breast cancers are detectable in histologically normal breast epithelium.

PURPOSE: Previously, we found that gene expression in histologically normal breast epithelium (NlEpi) from women at high breast cancer risk can resemble gene expression in NlEpi from cancer-containing breasts. Therefore, we hypothesized that gene expression characteristic of a cancer subtype might be seen in NlEpi of breasts containing that subtype. EXPERIMENTAL DESIGN: We examined gene expression in 46 cases of microdissected NlEpi from untreated women undergoing breast cancer surgery. From 30 age-matched cases [15 estrogen receptor (ER)+, 15 ER-] we used Affymetryix U133A arrays. From 16 independent cases (9 ER+, 7 ER-), we validated selected genes using quantitative real-time PCR (qPCR). We then compared gene expression between NlEpi and invasive breast cancer using four publicly available data sets. RESULTS: We identified 198 genes that are differentially expressed between NlEpi from breasts with ER+ (NlEpiER+) compared with ER- cancers (NlEpiER-). These include genes characteristic of ER+ and ER- cancers (e.g., ESR1, GATA3, and CX3CL1, FABP7). qPCR validated the microarray results in both the 30 original cases and the 16 independent cases. Gene expression in NlEpiER+ and NlEpiER- resembled gene expression in ER+ and ER- cancers, respectively: 25% to 53% of the genes or probes examined in four external data sets overlapped between NlEpi and the corresponding cancer subtype. CONCLUSIONS: Gene expression differs in NlEpi of breasts containing ER+ compared with ER- breast cancers. These differences echo differences in ER+ and ER- invasive cancers. NlEpi gene expression may help elucidate subtype-specific risk signatures, identify early genomic events in cancer development, and locate targets for prevention and therapy.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Cbf genes of the Fr-A2 allele are differentially regulated between long-term cold acclimated crown tissue of freeze-resistant and - susceptible, winter wheat mutant lines.

BACKGROUND: In order to identify genes that might confer and maintain freeze resistance of winter wheat, a comparative transcriptome analysis was performed between control and 4 wk cold-acclimated crown tissue of two winter wheat lines that differ in field freeze survival. The lines, generated by azide mutagenesis of the winter wheat cultivar 'Winoka' were designated FR (75% survival) and FS (30% survival). Using two winter lines for this comparative analysis removed the influence of differential expression of the vernalization genes and allowed our study to focus on Cbf genes located within the Fr-A2 allele independent of the effect of the closely mapped Vrn allele. RESULTS: Vernalization genes, (Vrn-A1, B1 and D1), and the transcription factor gene, TaVrt-2, were up-regulated to the same extent in FR and FS lines with cold acclimation thus confirming that azide mutagenesis had not modified the winter habitat of the lines. One category of Cbf genes, (Cbf-2, -A22 and B-22) reflected an increase in level of expression with cold acclimation in both FR and FS lines. Another category of Cbf genes (Cbf-3, -5, -6, -12, -14 and -19) were differentially expressed between cold-acclimated FR and FS lines relative to the non-acclimated controls. Comparison of expression patterns of the two categories of Cbf genes with the expression patterns of a set of ABA-dependent and -independent Cor/Lea genes revealed similar patterns of expression for this sample of Cor/Lea genes with that for Cbf-2 and -22. This pattern of expression was also exhibited by the Vrn genes. CONCLUSION: Some Cor/Lea genes may be co-regulated by the Vrn genes during cold acclimation and the Vrn genes may also control the expression of Cbf-2, -A22 and -B22. The increased freeze survival by the FR line and the increase in expression levels of wheat Cbf genes, Cbf-3, -5, -6, -12, -14 and -19 with cold acclimation in the FR line suggests a possible gain of function mutation resulting in higher levels of expression of these Cbf genes and increased freeze survival.

Gene expression profile analysis of human hepatocellular carcinoma using SAGE and LongSAGE.

BACKGROUND: Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide and the second cancer killer in China. The initiation and malignant transformation of cancer result from accumulation of genetic changes in the sequences or expression level of cancer-related genes. It is of particular importance to determine gene expression profiles of cancers on a global scale. SAGE and LongSAGE have been developed for this purpose. METHODS: We performed SAGE in normal liver and HCC samples as well as the liver cancer cell line HepG2. Meanwhile, the same HCC sample was simultaneously analyzed using LongSAGE. Computational analysis was carried out to identify differentially expressed genes between normal liver and HCC which were further validated by real-time quantitative RT-PCR. RESULTS: Approximately 50,000 tags were sequenced for each of the four libraries. Analysis of the technical replicates of HCC indicated that excluding the low abundance tags, the reproducibility of SAGE data is high (R = 0.97). Compared with the gene expression profile of normal liver, 224 genes related to biosynthesis, cell proliferation, signal transduction, cellular metabolism and transport were identified to be differentially expressed in HCC. Overexpression of some transcripts selected from SAGE data was validated by real-time quantitative RT-PCR. Interestingly, sarcoglycan-epsilon (SGCE) and paternally expressed gene (PEG10) which is a pair of close neighboring genes on chromosome 7q21, showed similar enhanced expression patterns in HCC, implicating that a common mechanism of deregulation may be shared by these two genes. CONCLUSION: Our study depicted the expression profile of HCC on a genome-wide scale without the restriction of annotation databases, and provided novel candidate genes that might be related to HCC.

Identifying nonspecific SAGE tags by context of gene expression.

Many serial analysis of gene expression (SAGE) tags can be matched to multiple genes, leading to difficulty in SAGE data interpretation and analysis. As only a subset of genes in the human genome are transcribed in a certain type of tissue/cell, we used microarray expression data from different tissue types to define contexts of gene expression and to annotate SAGE tags collected from the same or similar tissue sources. To predict the original transcript contributing a nonspecific SAGE tag collected from a particular tissue, we ranked the corresponding genes by their expression levels determined by microarray. We developed a tissue-specific SAGE tag annotation database based on microarray data collected from 73 normal human tissues and 18 cancer tissues and cell lines. The database can be queried online at: http://www.basic.northwestern.edu/SAGE/. The accuracy of this database was confirmed by experimental data.

Genome-wide analysis of antisense transcription with Affymetrix exon array.

BACKGROUND: A large number of natural antisense transcripts have been identified in human and mouse genomes. Study of their potential functions clearly requires cost-efficient method for expression analysis. RESULTS: Here we show that Affymetrix Exon arrays, which were designed to detect conventional transcripts in the sense orientation, can be used to monitor antisense expression across all exonic loci in mammalian genomes. Through modification of the cDNA synthesis protocol, we labeled single-strand cDNA in the reverse orientation as in the standard protocol, thus enabling the detection of antisense transcripts using the same array. Applying this technique to human Jurkat cells, we identified antisense transcription at 2,088 exonic loci of 1,516 UniGene clusters. Many of these antisense transcripts were not observed previously and some were validated by orientation-specific RT-PCR. CONCLUSION: Our results suggest that with a modified protocol Affymetrix human, mouse and rat Exon arrays can be used as a routine method for genome-wide analysis of antisense transcription in these genomes.